Dynamic characteristics of deep-sea ROV umbilical cables under complex sea conditions
Keywords
Nomenclature
Umbilical cable length
Number of discrete units
Umbilical cable diameter
Umbilical cable density
Sea water density
Umbilical cable cross-sectional area
Umbilical cable inertia tensor
Umbilical cable stiffness coefficient
Umbilical cable section position vector
Umbilical cable centerline arc coordinates
Force on the umbilical cable section
Moment of the umbilical cable section
- E
- G
T-axis unit vector
Umbilical cable section angle
Umbilical cable section velocity
Ocean current velocity
Tangential water resistance coefficient
Normal water resistance coefficient
Subnormal water resistance coefficient
Additional mass force coefficient
Umbilical cable gravity in water
Cage mass in water
Change rate of the section angular displacement relative to the arc coordinate
Section angular velocity
Distribution force
Force of the ocean current on the umbilical cable
Umbilical cable additional mass force
Umbilical cable gravity
1. Introduction
In recent years, abundant marine resources have gradually attracted the attention of all countries, especially with the depletion of resources on land, which has stimulated and driven the rapid development of deep-sea remotely operated vehicles (ROVs) (Chen and Liu, 2018). ROV systems are usually composed of a supporting mother ship, an umbilical cable, a tether-management-system (TMS, also named cage), and the ROV, as shown in Fig. 1. The umbilical cable is responsible for energy transmission, information interaction, and receiving and releasing the ROV, which is the lifeline of the ROV system (Li et al., 2013). Under complex sea conditions, the mother ship generates a large six-degrees-of-freedom motion in three-dimensional space. The heave motion of the mother ship has the most significant effect on the cage (Driscoll, 1999). The gravity and buoyancy of the ROV are almost equal in water, and its impact on the cage is negligible. The umbilical cable amplifies the motion amplitude of the mother ship owing to its elastic material, and a high cage motion amplitude may cause serious accidents, such as equipment damage and umbilical cable breakage. Therefore, theoretical research on umbilical cables is urgently required (Zhu et al., 2008). Fig. 1. Schematic diagram of a deep-sea ROV umbilical cable system.
Extensive research on the nonlinear mechanical model of ROV umbilical cables has been performed based on the mass-spring model, finite difference method and finite element method. Driscoll et al. (1999, 2000a, 2000b, 2000c) established a one-dimensional mass-spring model for umbilical cables, and the equation of motion was determined by assembling a finite element discrete force balance equation. Buckham et al. (2001, 2003) established a mass-spring model for an umbilical cable under low tension, and the internal bending force was deduced in terms of the local curvature. The mass-spring method is a low-order finite element method that has the advantages of modularity and can easily to embed cable modules with different properties. Park et al. (2003, 2005) used the finite difference method to model a towed sonar system, and considered the influence of the cable tangential and normal resistance coefficients. The numerical calculation results were in good agreement with the experimental data. Koh and Rong (2004) used the finite difference method and conducted a three-dimensional dynamic analysis of general cables used in engineering, considering the geometric nonlinear factors of stretching, bending, torsion and large deformation. The accuracy of the dynamic analysis was verified experimentally. Based on the finite element method, Cho et al. (2004) proposed a method for analyzing the dynamic response of an umbilical cable by considering the coupling between the mother ship and ocean current. Eidsvik and Schjolberg (2016, 2018) used the linear finite element theory based on the Euler–Bernoulli beam theory to construct a numerical model for an ROV umbilical cable, which is suitable for low-tension situations. Curic (2003) established an umbilical cable model with variable length based on the finite element method and evaluated its effectiveness using numerical simulations. Jordan and Bustamante (2007, 2008) ignored the gravity of the umbilical cable and analyzed the stability of the ROV system under the action of heaving motion and nonlinear vibrations. Quan et al. (2014, 2015, 2016, 2020), calculated the stress state of the ROV umbilical cable endpoint under the joint action of the mother ship and ocean current using the finite element method and the geometrically accurate beam model.
The elastic rod theory is a modeling method that is different from the mass-spring model, finite difference method and finite element method, and mainly includes the Cosserat theory (Liu and Xue, 2011; Kumar, 2016; Gao et al., 2017) and Kirchhoff rod theory (Kratchman et al., 2017; Luo et al., 2014; Bretl and Mccarthy, 2014). Compared with the Kirchhoff rod theory, the Cosserat theory considers the strain of the elastic rod, where the model is more accurate, but more complex and difficult to solve. Owing to its simple mathematical model, the Kirchhoff rod theory is widely used. Liu et al. (2018) established a curved surface constraint theoretical model for a soft cable based on the Kirchhoff rod theory. Liu et al. (2014) and Wang et al. (2012) used the Kirchhoff rod theory to develop a simulation for a cable in a virtual assembly process and obtained good simulation results. Based on the Kirchhoff rod theory, Goyal et al. (2005, 2008) studied the process of twisted submarine cables, and applied the model to DNA research. Currently, the application of the Kirchhoff rod theory in cable modeling mainly focuses on the simulation of the static form of the cable in the virtual assembly process. However, only a few applications in ROV umbilical cable analysis have considered the external environmental force and axial tensile deformation.
Based on the Kirchhoff rod theory, a nonlinear dynamic model is established, considering the axial tensile deformation, the force of the ocean current on the umbilical cable and the complex sea conditions. The remainder of this paper is organized as follows. Dynamic models of the umbilical cable are described in section 2. The static form of the umbilical cables with different lengths under different sea conditions and the dynamic characteristics under the joint action of the mother ship motion and ocean current are calculated in section 3. The conclusions are presented in section 4.
2. Dynamic analysis of an umbilical cable
2.1. Main model
A dynamic model of an umbilical cable can be built from the Kirchhoff rod theory in the world coordinate system , Frenet coordinate system , and spindle coordinate system , as shown in Fig. 1 (Jordan and Bustamante, 2007). The origin O of the world coordinate system is connected to the initial point of the umbilical cable centerline. The Frenet coordinate system follows the movement of point P. The axes are defined as follows: tangent vector , main normal vector , and vice normal vector . The spindle coordinate system is fixed to the cable section. The z-axis coincides with the T-axis, while the angle between the N and x axes is denoted by .
As shown in Fig. 2, the micro-arc is analyzed in the world coordinate system . The internal forces and torques in the negative section of point P are and , respectively. And the internal forces and torques in the section of point are and , respectively. The distributed forces, including the distributed forces of current and gravity, are f. At point P, we obtain(1)(2)where is the section rotation angle from time t to time . Fig. 2. Arc-microelement force diagram. The position of any point P on the centerline is determined by the vector r. The position of P and relative to point O are r and , respectively. The arc-coordinates relative to point O are s and , respectively.
From the momentum theorem and the moment of momentum theorem, we obtain(3)(4)where is the z-axis basis vector in the spindle coordinate system . The inertia tensor , and(5)
Let be the curvature-twisting of the section, which is the rate of change of the section angle relative to the arc coordinate s (Liu, 2006).(6)
and satisfy(7)
Taking the partial derivative of Eq. (1) with respect to s yields(8)
In the spindle coordinate system , Eqs. (3), (4), (7), (8) can be written as (Liu, 2006):(9)
Eq. (9) can be written in a scalar form as follows:(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20)(21)
The Euler angle is usually employed to describe the robot posture, which appears singular. The Euler parameters (,,,) can also be used to describe the posture of the section, which eliminates singular solutions. Therefore, the Euler parameters are adopted to describe the section posture of the umbilical cable. The relationship between the Euler parameters, and is defined as follows (Liu, 2006):(22)(23)
The partial differential equations for (,,,,,,,,,) can be obtained by unifying all the physical quantities with the Euler parameters, where the Euler parameters satisfy(24)
Before solving the dynamic model for the umbilical cable, the boundary conditions musted be defined. The initial point of the umbilical cable moves with the mother ship. Considering only the heaving motion of the mother ship and ignoring horizontal motion, the heaving motion of the mother ship at sea level can be expressed using a sine function:(25)where B is the amplitude and T is the period. The position of the umbilical cable at the initial point is(26)
The form of the section is limited by the Euler parameters. Supposing that the sections at the initial and end positions in the spindle coordinate system are set parallel to those in the world coordinate system , the Euler parameters at the initial and end position sections satisfy:(27)(28)
The force on the end section of umbilical cable is(29)
Eqs. (26), (27), (28), (29) constitute the boundary conditions of the umbilical cable dynamics.
Written in matrix form,(30)
The closed partial differential in Eq. (30) represents the umbilical cable dynamics model.
2.2. Distribution force
The rotation matrix of the spindle coordinate system with respect to the world coordinate system is expressed in terms of the Euler parameters as follows:(31)
The velocity of the current is approximately linearly proportional to the sea depth. There is almost no sea current below 2000 m. Therefore, it is assumed that the velocity is only a function of depth:(32)where is the current velocity at sea level and H = 2000 m. In the world coordinate system ,(33)
In the spindle coordinate system , the current velocity is(34)
In the spindle coordinate system , the relative speed of the ocean current and the umbilical cable is(35)
And the force of the ocean current on the cable is(36)
In addition, in a marine environment, the influence of the additional mass force on the umbilical cable model must be considered. The additional mass force of the cable in the spindle coordinate system is:(37)where is the density of sea water.
The umbilical cable is subject to gravity and a force distribution from the ocean current in the world coordinate system . Let the distribution of the gravitational force be . The projection of the force distribution in the spindle coordinate system is:(38)
2.3. Discretization of the equations
The differential quadrature method (DQM) algorithm is used to discretize Eq. (30) with the zero point of a Chebyshev polynomial as the node, as shown in Eqs. (39), (40).(39)(40)
Using the Lagrange interpolation basis function, the formula for calculating the weight coefficient matrix A is(41)where x represents an unknown quantity .
Owing to the tensile deformation of the umbilical cable centerline, the model must be modified when performing spatial dispersion. The arc coordinate in the stressed state is represented by , and the arc coordinate in the relaxed state is represented by s. The relationship between the two is given by(42)
The implicit time-domain algorithm is used to solve Eq. (31) in the time domain, as shown in Eqs. (43), (44) (Quan et al., 2014).(43)(44)(45)
2.4. Numerical solution
In Eq. (30), the partial differentials at each moment in time are combined to form a completely closed algebraic system of equations , where m is the number of equations. The problem can be transformed into a nonlinear least-squares problem, where the objective function is(46)
The optimization algorithm shown in Fig. 3 is used to solve the model. are preset parameters, k is the number of iterations, and is the initial value in each iteration. is the Jacobi matrix at , is the gradient, and GK is the modulus of .(47)(48)where is the search direction, which is divided into the Gauss–Newton (GN) and Levenberg–Marquardt (LM) directions. When the Hesse matrix is singular, a LM direction search is performed, and when the Hesse matrix is not singular, a GN direction search is performed. Sometimes, the determinant of the Hesse matrix is extremely large in the calculation process. Owing to the limited computing power of computers, the inversion fails. Therefore, when the Hesse matrix determinant is very large, the search is also performed in the GN direction.(49)(50) Fig. 3. Flow chart showing the optimization algorithm.
The Flag is an identifier, where indicates a search along the GN direction and indicates a search along the LM direction.
The optimal iteration step size is determined using the Armijo criteria.(51)
r is the ratio of the actual decline to the theoretical decline, and and are the thresholds of the trust-region radius. and are adjusted based on the search direction and trust region.(52)(53)
3. Calculation results
In this section, we verify the accuracy of the umbilical cable dynamic model using existing test data and then calculate the static and dynamic characteristics of the umbilical cable under different sea conditions. Data showing the motion of the mother ship and cage collected by the ROPOS Canadian deep-sea ROV system were used to verify the model (Driscoll et al., 1999). The relevant parameters of the umbilical cable are listed in Table 1. A comparison between the theoretical and actual values of the cage is shown in Fig. 4, for the case where the motion of the mother ship is considered and the effect of the ocean current is ignored. In Fig. 4(a), we calculate the movement of the cage based on the measured movement of the ROPOS's mother ship and compare it with the measured movement of the cage. As shown in Fig. 4(b), the theoretical calculation value of the cage motion coincides with the actual value, and the difference between the theoretical and the actual values was calculated. The mean difference is 0.0040. Table 1. ROPOS umbilical cable parameters. Fig. 4. Comparison between measured data and calculated values for the Cage. (c) means calculated values and (r) means real values. When the motion of the mother ship is ignored, the umbilical cable reaches static equilibrium under the action of the sea current, and total time derivatives in Eq. (30) are set to 0, that is, Eq. (30) is reduced to a steady state. Then, the variables in Eq. (30) become(54)(55)(56)(57)(58)(59)(60) Eqs. (54), (55), (56), (57), (58), (59), (60) constitute the static equations of the umbilical cable, which can be written in the matrix form as follows:(61) Before solving the static model for the umbilical cable, the boundary conditions must be defined. In the static model for an umbilical cable, the postures of the sections at the initial and final positions are the same as those in Eqs. (28), (29), respectively. However, Eqs. (26), (29) must be modified as follows:(62)(63) Eqs. (62), (63) define the boundary conditions of the static umbilical cable. The ocean current in the direction and there is no current below 2000 m. We calculate the static form of the umbilical cables under different lengths and sea conditions. When the umbilical cable is not affected by ocean currents, it remains upright. Under the action of ocean currents, the umbilical cable bends and deforms. We define the distance of the endpoint of the umbilical cable relative to the axis as the offset. The shape of the 2000 m long umbilical cable in ocean currents with different velocities is shown in Fig. 5. The umbilical cable has a larger deviation relative to the -axis as the ocean current gradually increases. When the ocean current velocities are 0.1, 0.2, 0.3, 0.4 and 0.5 m/s, the offsets of the umbilical cable end relative to the -axis are 1.326, 6.221, 13.367, 20.915 and 24.394 m, respectively, as shown in Table 2. At the same time, the equilibrium position of the endpoint of the umbilical cable gradually rises under the action of the ocean current, as shown in Fig. 5(c). However, the furthest point from the umbilical cable to the -axis is not the endpoint, but a point near the center of the cable. There are two reasons for this phenomenon: 1) the two ends of the umbilical cable are constrained, as shown in Eqs. (62), (63); 2) when the ocean current velocity decreases with depth, the force gradually decreases. Fig. 5. Umbilical cable shape at different ocean current velocities. Table 2. Offsets of the endpoints under different ocean current velocities. The shapes of the umbilical cables with different lengths when the ocean current velocity is 1.0 m/s are shown in Fig. 6. When the lengths of the umbilical cable are 1000, 2000, 3000, 4000, and 5000 m, the offsets of the endpoint to the z-axis are 14.412, 24.393, 17.828, 17.749 and 10.989 m, respectively, as shown in Table 3. The offset of the umbilical cable end first increases and then decreases because there is no current below 2000 m, and the umbilical cable only experiences gravity. Fig. 6. Shapes of the umbilical cables of different lengths. Table 3. Endpoint offsets for umbilical cables with different lengths. It is assumed that the direction of the ocean current is different in the range of 0–1000 m and 1000–2000 m. The velocity of the ocean current at sea level is 0.5 m/s, and the length of the umbilical cable is 2000 m. S1 ∼ S4 represent four cases where the directions of the current are different, and the direction of the current in each case is shown in Table 4. The shape of the umbilical cable for each case is shown in Fig. 7. When the direction of the ocean current is different at different depths, the static form of the umbilical cable is different, and the position of the endpoint is also different. The umbilical cable is bent in a three-dimensional space. Table 4. Direction of the ocean current. and represent the axes. Fig. 7. Shapes of umbilical cables with different current direction. The dynamic form of the umbilical cable under the joint action of the mother ship's motion and ocean current oscillates near the static location. The static shape can be used as a theoretical reference for the dynamic shape. At the same time, the form parameters can be used as the initial state of the algorithm in the dynamic analysis. Determining the dynamic characteristics of the umbilical cable under the joint action of the mother ship and ocean current is of great significance to the stability of the ROV system. This section describes the calculation of the dynamic characteristics of an umbilical cable under different conditions. Suppose the equation of motion for the mother ship is(64)where B = 1 m and T = 5 s. When L = 3000 m, we calculate the dynamic motion of the umbilical cables with different ocean current velocities. The static results are considered as the initial state of the umbilical cable, and we assume that each node has the same initial velocity but no initial acceleration.(65) As shown in Fig. 8(a)–(e), the phase difference gradually increases with an increase in the current velocity. When the ocean current velocities are 0.25 m/s, 0.5 m/s, 0.75 m/s, 1.00 m/s and 1.25 m/s, the corresponding amplitude amplification ratios are 1.530, 1.532, 1.538, 1.540 and 1.541, respectively, indicating that the amplitude amplification ratio gradually increases with an increase in the current velocity, as shown in Fig. 8(f). However, the force change at the head end of the umbilical cable is small. Compared to the umbilical cable gravity and cage gravity, the force generated by the ocean current is small, so the impact on the stress state of the umbilical cable is small. The umbilical cable has a stress of 80–160 kN at the head within a safe working load range. Fig. 8. Comparison of the motion and forces at the initial point with different ocean current velocities. In (a)–(e), the red dotted line represents the movement of the mother ship, and the solid blue line represents the movement of the end of the umbilical cable and the force at the head. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.) Suppose the motion amplitude of the mother ship is , and the ocean current velocity is . We calculate the effect of different periods of the mother ship on the umbilical cable. As shown in Fig. 9(a)–(e), with an increase in the motion period, the phase difference gradually decreases. When the mother ship motion periods are 2.5, 5.0, 7.5, 10.0, and 12.5 s, the corresponding amplitude amplification ratios are 1.56, 1.53, 1.51, 1.49, and 1.48, respectively. It can be observed that the amplitude amplification ratio tends to decrease with an increase in the motion period, as shown in Fig. 9(f). The stress on the umbilical cable decreases with an increase in the movement period of the mother ship. When , the force on the head of the umbilical cable is always within the safe working load range. However, when , the force on the head of the umbilical cable is greater than 200 kN, which exceeds the safe working load. Therefore, when the sea conditions are serious and the motion period of the mother ship is small, the umbilical cable is most likely to break. Fig. 9. Comparison of motion and force at initial point with different periods of mother ship motion. In (a)–(e), the red dotted line represents the movement of the mother ship, and the solid blue line represents the movement of the end of the umbilical cable and the force at the head. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.) Suppose the motion period of the mother ship is , and the current velocity is m/s. We analyze the effects of different mother ship amplitudes on the umbilical cable. As shown in Fig. 10(a)–(e), the phase difference gradually increases with an increase in the mother ship motion amplitude. When the mother ship motion amplitudes are 0.5, 0.75, 1.0, 1.25 and 1.5 m, the corresponding amplitude amplification ratios are 1.50, 1.52, 1.53, 1.54, and 1.55, respectively. The amplitude amplification ratio increases with an increase in the amplitude, as shown in Fig. 10(f). The stress on the umbilical cable increases with an increase of the amplitude. When the amplitude is 1.5 m, the maximum force on the head of the umbilical cable is 75–185 kN, which is within the safe working load range. Fig. 10. Comparison of motion and force at initial point with different amplitudes of mother ship motion. In (a)–(e), the red dotted line represents the movement of the mother ship, and the solid blue line represents the movement of the end of the umbilical cable and the force at the head. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.) Suppose the ocean current velocity is m/s, the motion period of the mother ship is , and the motion amplitude of the mother ship is . We calculate the dynamic motion of the umbilical cables of different lengths. As shown in Fig. 11(a)–(e), when the length of the umbilical cable increases, the phase difference increases. When the lengths are 1000, 2000, 3000, 4000 and 5000 m, the amplitude amplification ratios are 1.18, 1.32, 1.53, 1.47 and 1.43, respectively. As the length increases, the stress on the umbilical cable also increases. When the length is 5000 m, the maximum force on the head of the umbilical cable is 140–200 kN, which is within but very close to the safe working load range. Fig. 11. Comparison of motion and force at initial point when the mother ship is in sinusoidal motion with different cable lengths. In (a)–(e), the red dotted line represents the movement of the mother ship, and the solid blue line represents the movement of the end of the umbilical cable and the force at the head. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.) When the sinusoidal heave motion of the mother ship is transmitted to the cage through the cable, the heave motion of cage is also observed. Within a certain underwater depth range, the movement amplitude of the cage will be increased to the maximum, and the amplitude of the tension fluctuation of the cable will also be increased to the maximum. Determination of the depth range is essential for the safe operation of the deep-sea ROV systems. We call this depth range the resonance region of the deep-water umbilical cable system. As shown in Fig. 11(f), when the length is 3000 m, the amplitude amplification ratio is the largest, indicating a strong resonance phenomenon is generated at 3000 m. In this section, we use the motion data of the mother ship, which were measured by the "Da Yang Yi Hao" research ship. The measured data for the mother ship are shown in Fig. 12. The ocean current velocity is m/s. We calculate the dynamic motion of the umbilical cables of different lengths. Fig. 12. Real heave motion data of ship “Da Yang Yi Hao”. As shown in Fig. 13(a)–(e), when the length of the umbilical cable increases, the phase difference increases gradually. When the lengths of the umbilical cable are 1000, 2000, 3000, 4000, and 5000 m, the cage amplitude amplification ratios are 1.17, 1.26, 1.45, 1.43 and 1.39, respectively, as shown in Fig. 13(f). This trend indicates that the umbilical cable system exhibits a large resonance phenomenon at 3000 m. The heave amplitude of the cage is amplified underwater, but not infinitely. The cage amplitude amplification ratio varies from 1.1 to 1.5, in the range of 1000–5000 m. The stress fluctuation ranges at the initial point of the umbilical cable are 60–80, 80–110, 100–130, 130–160 and 150–190 kN, respectively. When the length is 5000 m, the stress of the umbilical cable is close to the safe working load of 200 kN. If the length of the umbilical cable continues to increase or the sea conditions become more adverse, the umbilical cable may break. A heave compensation device should be added to ensure the safe operation of the umbilical cable. Fig. 13. Comparison of motion and force at initial point when using mother ship actual motion data with different cable lengths. In (a)–(e), the red dotted line represents the movement of the mother ship, and the solid blue line represents the movement of the end of the umbilical cable and the force at the head. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.) We analyze the effects of different sea conditions, the movements of different mother ships and the length on the dynamic characteristics of the umbilical cables, as shown in Table 5. The results show that the amplitude amplification ratio is proportional to the current velocity and amplitude of the mother ship, and inversely proportional to the period of the mother ship. The force on the head of the umbilical cable is inversely proportional to the period of the mother ship motion and proportional to the length of the umbilical cable, and is less affected by the current velocity and amplitude of the mother ship motion. The phase difference is proportional to the current velocity, amplitude of the mother ship motion and length of the umbilical cable, and inversely proportional to the period of the mother ship motion. Table 5. Relationship between the dynamic characteristics of the umbilical cable and various parameters. CV: current velocity; PMM: period of the mother ship motion; AMM: amplitude of the mother ship motion; LUC: length of the umbilical cable; Based on the analysis presented in sections 3.3.4 Different umbilical cable lengths, 3.3.5 Actual motion action of the mother ship, it can be observed that when the length of the umbilical cable changes from 1000 to 3000 m, the amplitude amplification ratio of the cage gradually increases. When the length of the umbilical cable changes from 3000 to 5000 m, the amplification ratio of the cage amplitude gradually decreases, indicating that the umbilical cable system produces a strong resonance phenomenon in the 3000 m region. When it is far away from this resonance depth range, the heave amplitude of the cage is significantly reduced. We compare the effect of the different mother ship movements on the amplitude amplification ratios, as shown in Fig. 14. As the amplitude of the average motion of the mother ship increases, the amplitude amplification ratio also gradually increases. When the mother ship undergoes sinusoidal motion, the average motion amplitude and amplitude amplification ratio are the largest, which is the same conclusion as that in section 3.3.3. Fig. 14. Effects of different mother ship movements on amplitude amplification ratios.3.1. Model validation
symbol name value symbol name value d (mm) Diameter 30 Tangential water resistance coefficient 0.02 Gw (N/m) Cable gravity in water 25.9 Normal water resistance coefficient 2.0 mc(kg) Cage mass in water 4320 Subnormal water resistance coefficient 2.0 E (GPa) Young's modulus 64.4 Additional mass force coefficient 1.5 G (GPa) Shear modulus 26.9 Fs(kN) Safe working load 200 3.2. Static analysis of the mother ship's heaving motion with constant current
velocity 0.1 m/s 0.2 m/s 0.3 m/s 0.4 m/s 0.5 m/s offset 1.326 m 6.221 m 13.367 m 20.915 m 24.394 m length 1000 m 2000 m 3000 m 4000 m 5000 m offset 14.412 m 24.394 m 17.828 m 17.749 m 10.989 m Empty Cell S1 S2 S3 S4 0–250m 250–500m 3.3. Dynamic analysis under joint action
3.3.1. Different ocean current velocities
3.3.2. Different periods of the mother ship motion
3.3.3. Different amplitudes of the mother ship motion
3.3.4. Different umbilical cable lengths
3.3.5. Actual motion action of the mother ship
3.3.6. Discussion
Empty Cell CV↑ PMM↑ AMM↑ LUC↑ amplitude amplification ratio ↑ ↓ ↑ Take the maximum at 3000m force on the umbilical cable — ↓ ↑ ↑ Phase difference ↑ ↓ ↑ ↑
4. Conclusion
During the underwater operation of the deep-sea ROV, the umbilical cable often breaks owing to the joint action of the mother ship's motion and ocean current, resulting in the loss of the expensive underwater robot. Considering the shortcomings of existing modeling methods, a new modeling method based on the Kirchhoff rod theory was presented in this paper. After deriving a dynamic model, the DQM algorithm was used for spatial discretization, and the method was adopted for iterations in the time domain. The accuracy of the model was verified against the data collected by the Canadian ROPOS system. The static shape and the dynamic characteristics of the umbilical cables were calculated.
The motion amplitude and period of the mother ship and ocean current velocity have little effect on the amplitude amplification rate of the cage, but the length of the umbilical cable has a significant effect. In particular, when the length of the umbilical cable is 3000 m, the amplitude amplification rate reaches a maximum. Therefore, the cage should avoid staying in the 3000 m area for too long, as it may result in cable damage.
CRediT authorship contribution statement
Peng Chen: Modeling, Algorithm design, Writing – original draft, preparation. Yuwang Liu: Software, Data curation. Shangkui Yang: Visualization, Investigation. Jibiao Chen: Writing- Reviewing. Qifeng Zhang: Software, Validation. Yuangui Tang: Supervision.
Declaration of competing interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgments
This work was supported by the Key Research Program of Frontier Sciences (Grant No. ZDBS-LY- JSC011), National Natural Science Foundation of China (61821005, 51975566), National Key R&D Program of China (2017YFC030560101) and Liao Ning Revitalization Talents Program (XLYC1807090).
References
- Bretl and Mccarthy, 2014Quasi-static manipulation of a Kirchhoff elastic rod based on a geometric analysis of equilibrium configurationsInt. J. Robot Res., 33 (1) (2014), pp. 48-68
- Buckham et al., 2003Three dimensional dynamics simulation of slack tether motion in an ROV system2003 13th International Offshore and Polar Engineering Conference (2003), pp. 127-134
- Buckham and Nahon, 2001Formulation and validation of a lumped mass model for low-tension ROV tethersInt. J. Offshore Polar Eng., 11 (4) (2001), pp. 282-289
- Curic, 2003Development of a Variable Length Underwater Cable ModelMaster Dissertation, Florida Atlantic University (2003)
- Chen and Liu, 2018A novel voltage regulation strategy for the electric power delivery system of a 6000-m ROVAppl. Ocean Res., 80 (2018), pp. 112-117
- Cho et al., 2004Struct. Eng. Mech., 18 (4) (2004), pp. 429-440
- Driscoll, 1999Dynamics of a Vertically Tethered Marine PlatformDoctoral Dissertation, University of Victoria (1999)
- Driscoll et al., 2000aDevelopment and validation of a lumped-mass dynamics model of a deep-sea ROV systemAppl. Ocean Res., 22 (3) (2000), pp. 169-182
- Driscoll et al., 2000bThe motion of a deep-sea remotely operated vehicle system Part 1: motion observationsOcean. Eng., 27 (2000), pp. 29-56
- Driscoll et al., 2000cThe motion of a deep-sea remotely operated vehicle system Part 2: analytical modelOcean. Eng., 27 (2000), pp. 57-76
- Eidsvik and Schjolberg, 2016Time domain modeling of ROV umbilical using beam equations10th IFAC Conference on Control Applications in Marine System (2016), pp. 452-457
- Eidsvik and Schjolberg, 2018Finite element cable-model for Remotely Operated Vehicles (ROVs) by application of beam theoryOcean. Eng., 163 (2018), pp. 322-336
- Gao et al., 2017Mechanical model of dexterous continuum manipulators with compliant joints and tendon/external force interactionsIEEE-ASME Transactions on Mechatronics, 22 (1) (2017), pp. 465-475
- Goyal et al., 2005Nonlinear dynamics and loop formation in Kirchhoff rods with implications to the mechanics of DNA and cablesJ. Comput. Phys., 209 (1) (2005), pp. 371-389
- Goyal et al., 2008Non-linear dynamic intertwining of rods with self-contactInt. J. Non Lin. Mech., 43 (1) (2008), pp. 65-73
- Jordan and Bustamante, 2007Numerical stability analysis and control of umbilical-ROV systems in one-degree-of-freedom taut-slack conditionNonlinear Dynam., 49 (1–2) (2007), pp. 163-191
- Jordan and Bustamante, 2008Guidance of underwater vehicles with cable tug perturbations under fixed and adaptive control systemsIEEE J. Ocean. Eng., 33 (4) (2008), pp. 579-598
- Koh and Rong, 2004Dynamic analysis of large displacement cable motion with experimental verificationJ. Sound Vib., 272 (1) (2004), pp. 187-206
- Kratchman et al., 2017Guiding elastic rods with a robot-manipulated magnet for medical applicationsIEEE Trans. Robot., 33 (1) (2017), pp. 227-233
- Kumar et al., 2016A helical Cauchy-Born rule for special Cosserat rod modeling of Nano and continuum rodsJ. Elasticity, 124 (1) (2016), pp. 81-106
- Liu, 2006Nonlinear Mechanics of Thin Elastic Rod: Theoretical Basis of Mechanical Model of DNA(first ed.), Tsinghua University Press, Beijing (2006)(in Chinese)
- Liu et al., 2018Nonlinear mechanics of flexible cables in space robotic arms subject to complex physical environmentNonlinear Dynam., 94 (1) (2018), pp. 649-667
- Liu and Xue, 2011Stability analysis of helical rod based on exact Cosserat modelApplied Mathematics and Mechanics-English Edition, 32 (5) (2011), pp. 603-612
- Luo et al., 2014A robust and real-time vascular intervention simulation based on Kirchhoff elastic rodComput. Med. Imag. Graph., 38 (8) (2014), pp. 735-743
- Liu et al., 2014Physics-based modeling and simulation for motional cable harness designChin. J. Mech. Eng., 27 (5) (2014), pp. 1075-1082
- Li et al., 2013The tension calculation of self-propelled, cable-controlled vehicles' cablesAdv. Mater. Res., 823 (2013), pp. 24-28
- Park et al., 2003A numerical and experimental study on dynamics of a towed low tension cableAppl. Ocean Res., 25 (5) (2003), pp. 289-299
- Park et al., 2005A numerical analysis for the dynamic behavior of ROV launcher and umbilical cable under combined excitations2005 15th International Offshore and Polar Engineering Conference (2005), pp. 313-318
- Quan et al., 2016The nonlinear finite element modeling and performance analysis of the passive heave compensation system for the deep-sea tethered ROVsOcean. Eng., 127 (2016), pp. 246-257
- Quan et al., 2014Dynamics analysis of planar armored cable motion in deep-sea ROV systemJ. Cent. S. Univ., 21 (12) (2014), pp. 4456-4467
- Quan et al., 2015A geometrically exact formulation for three-dimensional numerical simulation of the umbilical cable in a deep-sea ROV systemChina Ocean Eng., 29 (2) (2015), pp. 223-240
- Quan et al., 2020Dynamics calculation for variable-length underwater cable with geometrically nonlinear motionOcean. Eng., 212 (2020), p. 107695
- Wang et al., 2012Dynamic simulation and disturbance torque analyzing of motional cable harness based on Kirchhoff rod modelJ. Mech. Eng., 25 (2) (2012), pp. 346-354
- Zhu et al., 2008A multi-body space-coupled motion simulation for a deep-sea tethered remotely operated vehicleJ. Hydrodyn., 20 (2) (2008), pp. 210-215
Cited by (14)
SPH-DEM modeling of cable-controlled ROVs: Underwater mobility and path planning
2024, Ocean EngineeringA hamiltonian global nodal position finite element method for dynamics analysis of submarine cables
2022, Ocean EngineeringCitation Excerpt :Yang et al. (2021) developed an analytical model to calculate the nonlinear tension-torsion coupling effect of the submarine cable. Chen et al. (2021) established the submarine cable dynamic model under complex sea conditions based on the Kirchhoff rod theory. Viel et al. (2022) proposed a model to predict the shape of the submarine cable so as to avoid entanglement of the cable itself or with surrounding obstacles.
Analysis of the vibration characteristics of non-uniform tension cables in a ship-cable-robot mining system
2023, Harbin Gongcheng Daxue Xuebao/Journal of Harbin Engineering University